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Abstract. We solve the problem of first-passage time for run-and-tumble particles in one dimension. Exact 
expression is derived for the mean hrst-passage time in the general case, considering external force-fields 
and chemotactic-fields, giving rise to space dependent swim-speed and tumble rate. Agreement between 
theoretical formulae and numerical simulations is obtained in the analyzed case studies - constant and sinu¬ 
soidal force fields, constant gradient chemotactic field. Reported findings can be useful to get insights into 
very different phenomena involving active particles, such as bacterial motion in external helds, intracellular 
transport, cell migration, animal foraging. 
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1 Introduction 

First-passage problems are of great interest in many disci¬ 
plines, from physics and chemistry to engineering and bi¬ 
ology [Dll]. After the seminal paper of Kramers on kinetic 
reactions [3] , a lot of efforts have been devoted to analyze 
first-passage processes, or the related narrow-escape prob¬ 
lems [1], in a variety of cases, both using theoretical mod¬ 
eling and numerical analysis [hllBllTlISlI^ ITOlITT] . The basic 
question is to find the mean first-passage time (MFPT), 
i.e. the mean time needed for a particles (that can be 
an atom, molecule, cell, animal or also a signal) to reach 
for the first time a specific site. Usually one considers 
processes in which the particle’s motion is described by 
Brownian dynamics and the probability density function 
(PDF) obeys a Fokker-Planck equation [T^]. This encom¬ 
passes many real processes involving thermal passive ob¬ 
jects, such as molecules, atoms or colloidal particles. In 
recent years a new kind of systems, characterized by an 
active nature, have been found to describe many interest¬ 
ing situations in real living world, from bacterial baths 
miM to cells migration m or animals’ movements CD- 
Such systems are characterized by self-propulsion mecha¬ 
nisms inducing a persistent motion, that can be described 
by simple random walk models [mill]. Among them the 
run-and-tumble model [14j has been found to be a very 
powerful modeling of real situations, allowing for analytic 
treatments [miiTiEi] and easy implementation in numer¬ 
ical simulations . In the simple version of the 
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model, particles perform a straight line motion at constant 
speed alternating tumble events (assumed to be random 
Poissonian processes), in which particle’s orientation ran¬ 
domly changes. The run-and-tumble model captures the 
behavior of motile bacteria, such as E.coli mmi, but is 
also related to many different physical processes, such as 
the motion of electrons in metals (Lorentz kinetic model) 
Eami the propagation of signals in transmission lines 
(telegraph equation) | 28 j . the dynamics instability of mi¬ 
crotubules mm. the Dirac equation for relativistic par¬ 
ticles in one spatial dimension mm- 
In this paper we analyze the first-passage problem for run- 
and-tumble particles in one dimension. We consider the 
general case in which external potentials and chemotactic 
fields are present, inducing a space dependence of parti¬ 
cles’ speed and tumble rate. A general expression of the 
mean first-passage time is given in the general case, and 
specialized to some interesting case studies - linear and 
sinusoidal potentials, constant gradient chemotactic field. 
Well known expression valid for Brownian particles is re¬ 
covered in the limit of high tumble rate and particle’s 
velocity. 


2 First passage time 


We consider run-and-tumble particles in one dimension. 
Denoting with P^{x,t) and P^{x,t) the PDFs of right- 
oriented and left-oriented particles, with and their 
swim speeds and with and the tumble rates, we can 













2 


L. Angelani, R. Di Leonardo, M. Paoluzzi: First-passage time of run-and-tumble particles 


write the continuity equations as follows 
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We treat here the very general case in which speeds (x) 
and tumble rates (x) are space dependent quantities, 
for example due to the presence of external force fields 
or chemotactic fields. We can write the equations for the 
total PDF P{x,t) = P^{x,t) + P^{x,t) and the current 
J{x,t) =Vj, {x)P^ {x, t) - {x)P^ (x, t) as 


The solution of Eq. Q can be written as 

2v,{L) 1 
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where, for the sake of simplicity, the dependence on space 
and time of the different quantities has not explicitly indi¬ 
cated. We want to determine the mean first passage time, 
that is the mean time needed for a particle starting its 
motion at x = 0 to reach for the first time the position 
L. In order to solve the problem we consider initial condi¬ 
tions P„(x, 0) = (5(x), P^ (x, 0) =0 and reflecting/absorbing 
boundary conditions respectively at x = 0 and x = L: 
J{0,t) = 0 and P^{L,t) — 0, i.e. J(P,f) =P(L,t)u^(P). By 
noting that the survival probability that the particle has 
not yet been absorbed at time t is F(t) = dxP(x,t) 
and the probability density of the first-passage time is 
(p(t) = —dtP, the MFPT r = J dt t (p{t) can be written as 


The mean first passage time r is finally obtained from Eq. 

©• . 

Some interesting limits are: 

i) The free particles case (absence of external force fields 
and chemotactic fields), i.e. v^ j^{x) = v and a„j^(x)=a, 
leads to 

where Dq = v^ja. It is worth noting that the MEPT 
can be written as the sum of a pure ballistic term plus a 
diffusive term. 

ii) In the case of no tumbling, ^^=0, we get 


'^no—tumble — 


dx 
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Prom Eq.s (^) the quantity Q(x) obeys the following 
equation 


with boundary conditions Q{L)v^^{L) = 1. It is worth not¬ 
ing that the above equation is similar to that obtained 
for the PDF of a run-and-tumble particle in the station¬ 
ary case [38l[2Tl[T^ . We have then reduced the problem of 
MFPT to the search of a solution of a stationary equa¬ 
tion with particular boundary conditions. We introduce 
the following quantities 


which is the time a particle with space-dependent velocity 
Vj^{x) takes to go form 0 to L. 

In the following we will consider the two cases in which 
particles are immersed in an external potential field (Sec¬ 
tion 3) or in a chemoattractant field (Section 4). A few 
case studies will be analyzed, deriving theoretical expres¬ 
sions and comparing them with numerical simulations. 
The latter are performed considering run-and-tumble par¬ 
ticles which obey the ID version of the Eq. of motion 
dtY = V e + y, i , with e the unit vector indicating the di¬ 
rection of motion of the particle, y the particles mobility 
and /(x) the external force. The direction e is updated 
with rate a — 7 e • Vc (c(x) is the chemotactic field) ran¬ 
domly choosing, in a uniform way, the new one - e = ±x 
in ID (see the following sections for details). 


3 External forces 

We first analyze the case in which particles feel an external 
potential field and chemotactic effects are absent, i.e. the 
tumbling rates of right-oriented and left-oriented particles 
can be considered equal and uniform 


.M) = ' 


(13) 


and the speeds are given by 


v^{x)=vP y /(x) 
v^{x)=v- y /(x) 


A{x) = exp 


(14) 
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with V the free-speed of the particles, fj, the particles mo¬ 
bility and f{x) = —dj;V{x) the forc e due t o an external 
potential V{x). In such a case, Eq.s become: 
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where D{x) is the local diffusion coefficient D{x) = [v^ — 

A^P{x)]la- 

It is worth noting that the Brownian case is obtained in 
the limit a,v —> oo with constant Bg = ja^ leading to 
the well known expression of the mean escape time for 
Brownian particles 


dx f dy (18) 

J X 

where (3 = y./Do = yia/v^. 

We consider here two kinds of external potentials. The 
first one is a linear potential 


1 

Do Jo 


V{x) = fox (19) 

where fg = A/L (A is the barrier height). In this case the 
mean escape time r has an explicit expression 


r = UfMo _ 1 


L 

t^fo 


( 20 ) 



where D = In the free limit, fo 0, one has 

T = Tfree " See Eq.( 111 - while for fg vjii one has —>■ 0 
and the MFPT diverges as t ~ {2la)eyi'p(yLID). In Fig. 
1 (upper panel) the MFPT r is reported as a function of 
/o for different values of L |37!. Theoretical predictions 
are in agreement with numerical simulations. 

In the Brownian limit, > oo at constant Dg = n^/a, 

one has 


Dg 


(pfoY 


DiJ.foL/Do _i\ _ ^21) 

^ ' M/o 


We observe that < r, indicating that in the Brownian 
limit particles take less time to reach absorbing bound¬ 
aries. It is worth noting, however, that the comparison is 
made considering equal diffusivities, and this does not im¬ 
ply that in real situation, at room temperature, thermally 
activated particles are faster than self-propelling ones. In¬ 
deed, just to make an example, considering parameters’ 
values suitable for E.coli (w~30 p,m s“^, p,~60 /rm pN“^ 
s“^, 1 s“^) one has that thermal diffusivity at room 

temperature is about four order of magnitude less than 
the self-propelling one D ~ 10^ yjo? s“^ , or, in other 
words, in order to have the same diffusivity one has to 
consider a temperature of 10® K. 


As second example we consider a sinusoidal potential 


V{x) 


A 

¥ 


(1 — costtx/L) 


( 22 ) 


Fig. 1. Mean first-passage time r of run-and-tumble particles 
in external force fields as a function of field strength. Upper 
panel refers to linear potential, Eq. (191, while lower panel to 
sinusoidal potential, Eq. (22l. Lines are theoretical predictions 
and symbols are from numerical simulations. Three different 
values of box length L are reported. In the upper panel dashed- 
lines correspond to the Brownian limit, Eq. (211. Quantities are 
expressed in reduced units 1321 . 


where A again is the barrier height. The corresponding 
force on particles is f{x) = —fg simrx/L , where fg = 
ttA/2L. We can write an explicit expression for the quan¬ 
tity A defined by Eq.(|^ 


A{x) = exp 


aL [W{x) - IT(0)] ] 




/ 


(23) 


where r] = y.fg/v and the quantity W is given by W{x) = 
tan“^ [r; cos( 7 ra;/L)/A/l — 77 ^]. The MFPT T can be ob¬ 
tained by numerical integration of Eq. (§, by using Eq.s 
( 9|10 ). In Fig. 1 (lower panel) the MFPT r is reported as 
a function of fg for different values of L m- It is worth 
noting that, as already stated above, one observes a diver¬ 
gence of the escape time r when the external force equals 
the self propelling one v/ji {fg —>■ 1 in reduced unit). In¬ 
deed, due to the finite velocity of run-and-tumble parti¬ 
cles, they can be easily confined by strong enough external 
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fields, inducing interesting trapping phenomena, recently 
investigated in random environments, such as speckle fields 
[55] or random distribution of obstacles [251131]. 


4 Chemotaxis 

We now analyze the effects of chemotactic fields. Bacte¬ 
ria, such as E. coli, are able to sense the environment and, 
through a chemical internal circuit, modify their tumble 
rate in order to achieve a net preferential motion towards 
high nutrient concentration m- The adopted strategy 
consists in evaluating changes in ambient chemical con¬ 
centration performing a time-integral of concentration c 
filtered by a suitable kernel function [13] . A simplified ex¬ 
pression of the tumble rate a as a function of chemotac¬ 
tic field c(x) is obtained in the limit of weak concentra¬ 
tion gradient, giving rise to space dependent right and left 
tumble rates 

a„(x) = a — jv dxc(x) , , 

^ / (24) 

(x) = a + jv OxCiyX) 



Fig. 2. Mean first-passage time r of run-and-tumble parti¬ 
cles in chemotactic fields as a function of field strength. The 
two curves refer to the cases of chemoattractant boundary, 
Co > 0 (full symbols) and chemorepellent boundary, co < 0 
(open symbols). Lines are theoretical predictions and symbols 
are from numerical simulations. Quantities are expressed in 
reduced units I3Z|. 


where a is the tumble rate in the absence of chemoattrac¬ 
tants and 7 the strength of the particles reaction to the 
chemicals [Il|33|. It is worth noting, however, that mod¬ 
eling chemotaxis through spatially varying tumble rates 
could be a quite crude approximation of the mechanism 
adopted by real organisms, which instead involve more 
complex time integration m- By considering force-free 
particles 

Vn.L{x)=v (25) 

we have the following expression for the quantities A, B, 
and Q 


A{x) 

B{x) 


exp{/ dy-f dyc{y) 


dy a [vA{y)\ 
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and the MFPT takes the form 


(26) 

(27) 

(28) 


is at the absorbing boundary L, resulting in a net pref¬ 
erential swimming towards such a point. For cq < 0 the 
maximum is at x = 0 and reaching the boundary is made 
more difficult by the chemotactic field. For weak concen¬ 
tration gradient one has: 


T _ Tfree Cq ( ^ + 3^ ) 


(32) 


which, in the absence of chemicals, cq = 0, reduce to the 
free case t = T/ree- In Fig- 2 the MFPT r, Eq. ([m]), is 
reported as a function of |co| for the two cases of positive 
and negative cq, corresponding, respectively, to chemoat¬ 
tractant and chemorepellent boundary. 

We conclude this section discussing a different form of 
chemotaxis, in which the cell sets its tumble rate consid¬ 
ering the absolute attractant concentration instead of its 
temporal changes. This is the case, for example, of mutant 
bacteria lacking particular enzymes required to perform 
temporal comparisons of concentration [33] . If we assume 
instantaneous sensing we can write tumble rates as 


f dx 

pL 

^--yciL) a / g-7c(!/) 

Jo V 

V Jx 


(29) 


0!n,L(x) = a-7 c(x) 

leading to the following expression for MFPT 


(33) 


We consider here the simple case of constant gradient con¬ 
centration: 

c(x) = Co X (30) 


The MFPT reads 
L 


1 


(1_ !(I^)(l_e-T"=o^) (31) 

a 


ICqDq 

For Co > 0 the maximum chemoattractant concentration 


Tfree - V / dx X c(x) 


(34) 


which, in the case of linear concentration c(x) = co x, 
reduces to 


^ - yree 3^3 


(35) 


However, real cells measure chemical concentrations in a 
finite time, and then a more appropriate form of tumble 
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rates in the case of mutant bacteria is [33] 

a^(x) = a — 7c(x) -I- dxc{x) 

a^(x) = a — ^c{x) — jv dxc{x) 


(36) 


The ratio 7/7 is proportional to the average time the cell 
spends to measure chemical concentration. The MFTP 
now reads 


r = 


dx 


,—yc{x) 


o7c(i) 


dy 
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g7c(y)(Q, _ ^c{y)) 


(37) 


and, for constant gradient field c(x) = cq x 
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For weak concentration gradient one has 


(38) 


— '^free Cq 


7 ^ 

2v 


l + 7(a-2i) 

6V 7 


(39) 


indicating a reverse taxis (r > Tfree) when 7/7 < 3v/2L+ 
a/2, i.e. when the average time required for the cell to 
measure chemical concentration is long enough [33] . 


5 Conclusions 

First-passage time problems for active particles ruled by 
run-and-tumble dynamics are solved in one dimension, 
considering generic space-dependent external forces and 
chemotactic fields. Exact expressions are given for the 
mean first-passage time and tested with simulations in 
simple case studies. A diffusion limit is obtained for high 
particle’s velocity and tumble rate at constant diffusiv- 
ity, where well known expressions for Brownian particles 
are recovered. Our results can be of interest to many dis¬ 
ciplines, allowing a better understanding of phenomena 
involving persistent motion, such as barrier escaping of 
active particles intracellular transport forag¬ 

ing and predation phenomena in the animal kingdom [431 
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